Detection of single electron heat transfer statistics 
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We consider a quantum dot system whose charge fluctuations are monitored by a quantum point 
contact allowing for the detection of both charge and transferred heat statistics. Our system consists 
of two nearby conductors that exchange energy via Coulomb interaction. In interfaces consisting 
of capacitively coupled quantum dots, energy transfer is discrete. In the presence of hot spots, we 
predict gate dependent deviations away from a fluctuation theorem that holds for charge currents 
in the absence of temperature gradients. A fluctuation theorem holds for coupled dot configurations 
with heat exchange and no net particle flow. 
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I. INTRODUCTION 

Contrary to charge, heat currents are difficult to mea- 
sure on nanoscopic scales^. Seebeck-Onsager coefficients 
relate heat flows to voltage drops^, allowing for the detec- 
tion of mesoscopic heat currents in the linear regime"^'^. 
Recently multiterminal arrangements have been used to 
inject current and as thermometers in heat transistors^. 
However, a method for a direct measurement of non lin- 
ear electronic heat currents and their fluctuations has not 
been described. It would allow for instance to investigate 
recent proposals of refrigerators^, efficient heat convert- 
ers^"—, rectifiersi^ or diodes^ based on single electron 
transistors. 

Measurements of the charge counting statistics have 
become a standard way to access the noise properties 
of the electric flows through quantum dot systemai^"— . 
There, time resolved charge fluctuations are detected by 
current changes though a capacitively coupled quantum 
point contact (QPC). Experimental precision makes pos- 
sible to measure the behaviour of high order cumulants 
or finitie ferquency statisticaiiii^. For this reason, meso- 
scopic conductors offer an ideal playground for experi- 
mental verifications of nonlinear fluctuation relations as 
expressed for electric conductioiJ^"— . Pioneering experi- 
ments have been carried out recentlj*^^^. In experiments 
where the detector is asymmetrically coupled to a dou- 
ble quantum do t^^i^^ , different charge distributions are 
resolved as different outputs in the QPC. This way tra- 
jectories of charges flowing along and against an applied 
voltage can be traced out. If the external force is thermal, 
so temperature is inhomogeneous, the fluctuation theo- 
rem for charge currents does not hold. One must then 
consider both charge and energy currents^'' or, alterna- 
tively, heat. To test fluctuation relations, the ability to 
measure heat currents is therefore also of fundamental 
importance. 

Here we show that charge counting statistics measure- 
ments can give access to heat flows and their fluctuations 
in systems of quantum dots. It allows for tests of non lin- 
ear fluctuation relations for transport in the presence of 
thermal and voltage gradients. To illustrate this gener- 



ally, we consider a model of capacitively coupled conduc- 
tors that allows to separate heat and charge flows: charge 
flows through the terminals of each conductor; energy is 
exchanged via electron-electron interactions, as sketched 
in Fig. [TJ Two quantum dots are close enough such that 
their charge is mutually affected by Coulomb interac- 
tion with a capacitance coupling, C. Experiments have 
demonstrated a strong interdot coupling with no electron 
transfer through the capacitor— i^. Each quantum dot 
is coupled by tunneling barriers to different electronic 
reservoirs considered at local equilibrium. One of them, 
that we denote as s, allows charge transport between two 
terminals at the same temperature but different voltages. 
The other one, g, which is coupled only to one reservoir 
at a different temperature, supports no stationary charge 
but heat current. In previous work, we have shown that 
such geometry is able to generate a charge current in 
the unbiased two terminal conductor by converting the 
heat flowing though the gate quantum dolj^i^S. Energy 
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FIG. 1: Proposed scheme for heat detection. Two capacitively 
coupled quantum dots are connected to three terminals. One 
quantum dot (s) supports charge transport between two ter- 
minals. The other one {g) leads to heat transfer from a hot 
reservoir. A quantum point contact is asymmetrically coupled 
to both dots serving as a four state charge detector. Time re- 
solved cycles involving all the states, 1— >-2— >-3— ^4, represent 
the transfer of one quantum of heat, Ec, from the hot to the 
cold conductor. The reversed sequence, 1'— i>2'— >-3'— >4', takes 
heat from the cold to the hot reservoir. 
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dependent tunneling barriers provide the necessary sym- 
metry breaking"^^. The perpendicular directions of heat 
and charge currents is an important issue for energy har- 
vesting from an environment. Related effects have been 
predicted in a variety of configurations, when the third 
hot terminal consists of a phonon batb^>^, the electro- 
magnetic environment^, incoherent radiation™ or a spin 
insulator-ii. Recently, current generation from absorp- 
tion of external radiation has been detected in periodic 
semiconductor heterostructures^^*^, following ideas by 
one of us^^i^ and van Kampen^. 

Placing a QPC close to the quantum dots allows to de- 
termine the charge occupation of the system. The trans- 
parency of the constriction is affected by the electrostatic 
coupling to the charge of the quantum dots. In general, 
the coupling of the QPC to each quantum dot will be dif- 
ferent, strongly depending of their distance. This makes 
the current flowing through the QPC sensitive to the dif- 
ferent charging states. We emphasize that similar ar- 
rangements have allowed directional charge countin g^^i^^ 
or the detection of quantum operations^' in double quan- 
tum dot systems. Alternatively, more complex setups of 
two QPCs, each one coupled to one quantum dot, can 
be used^. The backaction of the detector on the system 
charge states can be tuned to be negligible. 

We propose detection of energy flows based on state re- 
solved charge counting. We consider a QPC detector with 
no backaction on the system. In our case, we will focus on 
a Coulomb blockade regime where each quantum dot is 
described by a single discrete level that can be occupied 
by up to one electron. That reduces our configuration 
to a four state basis: {Ns, Ng)^{{Ofl), (1,0), (0,1), (1,1)}. 
The generalization to cases with several levels is straight- 
forward. Up to sequential tunneling, three states only are 
not enough to induce correlations between the two con- 
ductors^. Due to level discretization, interdot energy 
flows occur in quanta, Ec=2q^/C, defined in terms of 
the geometrical capacitance of the coupled system, G^. 
Sequences can be identified in the time trace of the cur- 
rent through the QPC, (0,0)^(1,0)^(1,1)^(0,1)^(0,0), 
that represent the transfer of an energy Ec from the hot 
gate to the conductor, as sketched in Fig. [T] The re- 
versed trajectory transfers energy from the cold to the 
hot system. 



II. MODEL 

We are interested in a system supporting non parallel 
electronic charge and heat currents. The three termi- 
nal device described above is the minimal model of this 
characteristics. Heat is transferred between the two dots 
by electron-electron interaction only. At sufficiently low 
temperatures, phononic heat transport can be neglected. 
Coulomb interaction is modeled by the capacitances as- 
sociated to the tunneling barriers, Ci and the interdot 
coupling, C. In a self-consistent treatment, they de- 
fine the charging energies, Ua.ni{Vi}), of quantum dot 



a when the occupation of the other dot is n = 0, 1^. The 
quantum of transferred energy is given by the interdot 
charging energy, Ec—Ua,i-'Ua,o—'^'}^/C, with the effec- 
tive capacitance C={C-ssCj:g—C'^)/C, C^s=Ci+C2+C, 
and Csg=Cg+C. The three terminals are assumed to be 
at thermal equilibrium with voltage Vi and temperature 
T/. We will focus on the case that the two terminals of 
the conductor are at the same temperature, Ti=T2=Ts 
but can be at different voltages. 

A. Master equation 

The dynamics of our system can be described by a mas- 
ter equation for the reduced density matrix of the quan- 
tum dot system, p. We consider the sequential regime, 
where only diagonal elements of the density matrix are 
important^i^. Processes in which an electron tunnels 
out (-I-) or into (— ) a quantum dot through junction I 
with 71 = 0, 1 electrons in the other quantum dot are de- 
scribed by the tunneling rates T^. They are calculated 
by Fermi's golden rule, giving: T~^rinf[{Ean-qVi)/3i], 

rto=r;n-r;~, with Pl = {kTiy^, Ean^Sa + Uan, being 

f{x) = {l+e^)~^ the Fermi function. Sa is the bare en- 
ergy of the discrete level in quantum dot a. The rates 
Tin are in general energy dependent^^— . In particular, 
they depend on the charge occupation of the coupled dot. 
In the regime kT^Tin, broadening of the energy lev- 
els can be neglected. Written in matrix form, such that 
P—{poQi Pw, P017 Pii)j the master equation reads p=M.p, 
with 
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and T^=T^+T^. For convenience, we can sep- 
arate into diagonal and non diagonal terms: 
M.=AiD+ J2i n k ^in- matriccs J'l^ correspond to 
those jumps described by tunneling rates F^^. They al- 
low one to express the charge, energy and heat stationary 
currents, /c,i, /e,i and /h,i, respectively, in a compact 
form: 

(/a,z)=^tr[er„(j-+- (2) 

n 

where represents the transported amount in each 

tunnefing event: Qf^—I i^^"^ charge), 9^=i?Q„ (for en- 
ergy) and Qf^-^=Ean~qVi (for heat). Note that they are 
related by: /h,i=^e,i— ^c,i^/- We obtain p from the 
steady state solution given by M.p — 0. From charge 
conservation, we know that /c,g=0 and Ic,s—Ic,2=~Ic,ii 
which allows us to define the total charge current flowing 
through the conductor, Iq^s- The state resolved charged 
currents, (Iin), are given by the terms inside the sum in 
Eq. @. 
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III. COUNTING STATISTICS AND 
FLUCTUATION THEOREMS 

Briefly, in the master equation formalism the cumulant 
generating function of the counting statistics is obtained 
from the eigenvalue closest to zero of the matrix Ad-^, as 
determined from Al when introducing the counting fields 
Xa,i multiplying the corresponding off diagonal elements, 
j;^, with z^;,^=e*''®°"''".'i^. Different counting fields, 
Xa, with a = C,E, H, give rise to the statistics of the 
different charge, energy and heat currents. Doing so, 
and noticing that F^^ = e^'^'^Tj^^, which translates to 

J,t = e^^^^^{J,-f, (3) 

one can easily verify the symmetries of the characteristic 
polynomial det| Al,^;— A1|=0. For uniform temperatures, 
Ti—T,yi, one obtains the fluctuation theorem for charge 
transport: 

J'iiXc)=J'i-iXc + Ac), (4) 

with the affinities Ac i=Vif3r^. At zero applied voltage, 
one expects that detailed balance restores equilibrium. 
However, a finite charge current flows in the presence of 
hot spotsfi.!^. This apparent broken detailed balance is 
restored when considering the energy flows as well by in- 
troducing the energy counting fields. For that one has to 
include the product Zq j^^Zg into the master equation. 
By doing so, the combined charge and energy cumulant 
generating function is invariant under the transformation 
{^Xc^Xe} ->■ {-iXc+Ac, -iXE-Au}, with the en- 
ergy affinities AH,i=/3f^. Notice that the local detailed 
balance condition expressed in Eq. (|3]) involves the heat 
carried by electrons in the tunneling processes. Thus a 
more compact form of the fluctuation theorem can be 
written in terms of heat currents — which we showed to 
be related to charge and energy currents — when consid- 
ering the heat counting fields j : 

J^iiXH)=J'{-iXH-AH), (5) 

valid in the presence of both voltage and temperature 
gradients. As mentioned above, a way to gather informa- 
tion about the heat transfer is by state resolved count- 
ing statistics, once the energies Of„ are known. In our 
configuration, this information can be accessed experi- 
mentally. Furthermore, if we are interested in the heat 
transferred between the two conductors, only the differ- 
ence Qf^—QfQ—Ec becomes relevant. 

A. State resolved fluctuation relations 

If one is able to distinguish what state of the quan- 
tum dot system is involved in the tunneling of the de- 
tected charges, the same information is given by state 
resolved counting. This is the same as introducing 
zj^^= exp[ikQn] multiplying the corresponding rates F^^ 



in the off-diagonal elements of the matrix Ai. Nin de- 
notes the number of transferred charges through terminal 
I and from state n, and Qn are the level resolved parti- 
cle counting fields. Note that Ni—J^n-^in- The explicit 
level dependence of local deltailed balance conditions in 
Eq. ([3]) provides the fulfillment of the state resolved fluc- 
tuation relation: 

T{tC)=T{-tC-A), (6) 

where C and A;^„=0[^„/3; are the state resolved charge 
counting fields and affinities, respectively. The dimen- 
sion of the vectors C and A is the number of termi- 
nals times the number of states. In contrast to Eq. (jH), 
the relation shown in Eq. (jG)) is valid also in the pres- 
ence of temperature gradients. As a drawback, it is 
no longer universal but depends on the internal con- 
figuration of the system — in particular, on gate volt- 
ages. However, we stress that, in a self-consistent treat- 
ment, Qf„=Ea^n{{Vm})—qVi depcuds only on differences 
of voltages, so Eq. ^ is gauge invariant. 

In a recent work, Krause et alM^ showed that state 
resolved counting of electrons allows ffuctuation relations 
to be written when the statistics of energy dissipated is 
also measurable. That is not the case for phonons. In our 
setup, the energy absorbed by a conductor is only due to 
electron-electron interactions and can be treated in the 
same footing as the particle currents, as we show below. 
State resolved counting of electrons therefore provides 
the required information about the energy transfers. 

Our system allows for such state resolved counting. 
The quantity of interest is the heat flowing through the 
gate. This is independent of the number of levels of each 
dot. The relevant index is given by the charge occupa- 
tion. For quantum dots that permit up to one electron, 
n—1, labels the charge occupation of the dot that is not 
involved in the tunneling event. The conductor absorbs 
a quantum of energy Ec from the gate when an electron 
enters the gate dot occupying the state (1,1) and leaves it 
from state (0,1) (after the electron occupying the conduc- 
tor dot has been transferred carrying the extra energy to 
any lead). Thus, the amount of energy transferred to the 
gate lead is proportional to the difference of the number 
of electrons that have tunneled from the two states (1,1) 
and (0,1): 

Eg=Ec{Ng,-Ngo). (7) 

This equation relates the transferred heat statistics to the 
charge fluctuations in the gate, which can be measured 
in a state dependent detector as the one shown in Fig. [T] 
Charge and energy conservation can be written in 
terms of state resolved counting: ^jg^A^/=A^g=0 and 
J2ies(.^i;i~^i;o)=^g,i-^g,o, respectively. They pro- 
vide additional symmetries of the cumulant generating 
function expressed as shifts of the counting fields. If we 
separate the contributions from each conductor, system 
and gate, C—{Csi Cg)i we obtain: 

J^{iCs,iCg)^J'{iCs + a,iCg + b), (8) 



4 



for arbitrary parameters a and h. If, on the other 
hand, we separate particles tunneling at different ener- 
gies, = (Oo,Oi): 

^{KsO, iCsl,iCgO, <gl) = 

T{iCsO-a, iCsi-b, iCgO+b, <gi+a)- (9) 

The dimension of vectors (^sn and Cgn is given by the 
number of levels in each conductor times the number of 
terminals to which it is connected. Eqs. ^ and © lead 
to an incomplete fluctuation theorem when only electrons 
flowing through the conductor are measured^: 

T{iCsn)=J^i-iCsn-Asn-Agfi). (10) 

Note that the system and gate affinities are taken at dif- 
ferent occupations, n=not(n). It can be useful in cases 
where the QPC is not able to distinguish the states (0,0) 
and (0,1). 

The state resolved cumulant generating function 
is defined in terms of the probability to detect N 
transferred electrons during a period of time, t: 
J^(i^)= In P(7V, t)e~'''^^. For later convenience, the 
charge currents will be written in terms of the number of 
electrons transferred in a given time: Iin=qNin/t. Then, 
in the long time limit, Eq. (jlO[) can be expressed in terms 
of the probability to measure a given current configura- 
tion. Ism in the conductor—: 

l£s,n 

In the case of interest here where the two terminals of 
the conductor s are at the same temperature, the pre- 
vious expression can be nicely rewritten in terms of the 
total charge current, /c,s='? X^n ^2n/i, and the amount 
of energy absorbed from the gate, lB_,g—Eg/t, expressed 
in terms of particle flows in the conductor, cf. Eq. ([7]): 

^ln;^^-/c,s(l^i-V2)/?.-/H,g(/3,-/?.)/2. (12) 

It resembles fluctuation theorems for thermoelec- 
tric transport that separate thermal and electric 
forces^Li^Ti^. However, Eq. is configuration de- 

pendent: the right hand side depends on state resolved 
currents through /n.g- For instance, it is affected by 
gate voltages. For a single level quantum dot, it reads: 
lH.g=2Ec{Ic,s—Iii+l2Q), where the last two terms rep- 
resent the difference of state resolved currents /;„ flow- 
ing through the different leads (1=1,2) of the conductor 
at the different charge occupations (n=l,0) of the gate 
quantum dot. 

An illustration of the previous expression is plotted 
in Fig. [5] We denote the right hand side of Eq. 
with ^. Remember that Jn.g can be written in terms of 
charge currents. Remarkably, in asymmetric configura- 
tions that manifest noise induced transport, (^) does not 
vanish even if the total current, (/c.s), does at a certain 
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FIG. 2: Fluctuation relation. We plot the right hand 
side of Eq. H12p written in terms of mean charge cur- 
rents: (e) = (/c>(K-V2)/3.-Sc((/c,=>-(/ii) + (/2o))(/3<,-A), 
as a function of the bias voltage applied to the con- 
ductor. The inset shows the same curves normalized 
to the total charge current, (/c,s), for clarity. Differ- 
ent tunneling rate configurations are shown for homoge- 
neous {Ts=Tg) and inhomogeneous (Tsj^Tg) temperatures: 
"symmetric" (Fjn^r, V{?,n}), "asymmetric" (ri„=r, except 
rii=r2o=r/10), and "optimal" (r,„=r, except rii=r2o=0). 
The temperature of the conductor is kept kTa=5hT, while the 
gate is heatted to kTg=10hr in the inhomogeneous case. The 
presence of a hot spot generates a shift of the minimum. For 
asymmetric configurations, it does not cross the origin ex- 
cept for the optimal case, in which it goes to zero at the stall 
potential of Eq. (fT5)) . 

stall potential. Though at that voltage no net charge 
currents flow, heat flowing through the gate drives the 
system out of equilibrium. Detailed balance is restored 
only at the optimal converter configuration, cf. Eq. (fT3)) 
below. 

Optimal heat to charge conversion occurs in particu- 
lar cases where the tunneling events (0,0) -H- (1,0) and 
(0,1) ^ (1,1) can only take place through different ter- 
minals of the conductor—. For instance, if Fn-cFio and 
r2o^F2i, so 111, /2o^'0. Then, a universal relation holds: 

^'^]?p^=^c,s[(V^i-1^2)/3s-i?c(/3g-/3.)]. (13) 

The involved cycles are illustrated in Fig. [T] At this 
configuration, a charge q is transported across the con- 
ductor for each quantum of heat Ec that is exchanged 
between the two systems, so Ic,a/q—lYL,g/Ec- This im- 
plies that /c,s and /u.g are maximally correlated, i.e. 
{Ic,sI'H,g) = {Ic,s){I'H.g), and all their cumulants coincide. 

Note that a related expression can be obtained using 
the relation of rates for the cycle that transfers an energy 
Ec to the conductor, 7~'"=Fj~QF;^]^rJ^rjQ, and its reversed 
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FIG. 3: Heat diode in a two terminal configuration. Here, 
no charge but only heat transport takes place through the 
capacitor. 



sequence, 7 ^rgpFj^Fi^j^rj^Q, which gives: 
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■=exp[g(yi-F2)/3.+£;c(/3.-/3ff)]. 



(14) 



The exponent in the right hand side of Eq. (jl4p is given by 
the entropy flux per cycle^i^. We can therefore interpret 
the right hand side of Eq. (fT3|) : it relates the entropy flux 
to the Joule heat in the conductor and heat absorbed 
from the gate. It also provides the condition for the stall 
potential 



Ps 



(15) 



where a conversion efficiency approaching the Carnot 
limit is predicted^. 



IV. PURE HEAT TRANSPORT 

As discussed above, charge current statistics suffer de- 
viations from a universal behaviour in the presence of 
thermal gradients — unless one has access to the whole 
level dependent transferences, distinguishing the involved 
terminal as well. On the other hand, an applied voltage 
to the conductor prevents a fluctuation relation in terms 
of the gate charge only. One can avoid such difficulties 
by considering a two terminal system where no stationary 
charge current will flow^^. A schematic representation of 
such a system with the QPC detector added is shown in 
Fig. [3] Note that in this configuration, heat and energy 
currents coincide, lH,i=lE.\—EcIii/q, and therefore are 
conserved: Z); -^h,i= Z); -^e,i=0. 



In this case, the validity of the fluctuation theorem 
can be verified by detecting the statistics of charge trans- 
ferred from a given energy. The constraints introduced by 
charge conservation in each conductor {N1Q——N11) and 
energy conservation {Nsi=—NgQ) permits to write sep- 
arate relations for particles emitted through any state 
or terminal, Nin- For example, one can simply count 
the number of electrons Nsn emitted to the cold system 
conditioned on the hot dot containing n electrons. The 
corresponding cumulant generating function obeys the 
symmetry: 



(16) 



For particles emitted through the gate terminal, replace 
s ^ g in the previous expression. No gate voltage de- 
pendence enters Eq. dH 



V. CONCLUSIONS 

To summarize, we propose a mechanism to detect the 
statistics of electronic heat currents in Coulomb block- 
aded quantum dot systems. Multilevel charge detection 
by means of a side coupled quantum point contact al- 
lows to count the transfer of discrete energy Ec- The 
same scheme permits to investigate deviations of charge 
flow statistics in the presence of hot spots, manifested 
in non universal fluctuation relations. We investigate a 
three terminal heat to charge converter where noise in- 
duced transport avoids detailed balance of charge flows. 
Only for an optimal configuration when the absorption 
of a quantum of heat Eq leads to the transfer of a charge 
q across the conductor, can a fluctuation theorem for 
total charge flow be written. In a two terminal device 
with no net charge currents, state resolved charge count- 
ing statistics fulfill the fluctuation theorem as obtained 
for heat currents. All our results can be probed within 
present experimental technology. 
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